Vibrational modes of metal nanoshells and bimetallic core-shell nanoparticles 
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We study theoretically spectrum of radial vibrational modes in composite metal nanostructures 
such as bimetallic core-shell particles and metal nanoshells with dielectric core in an environment. 
We calculate frequencies and damping rates of fundamental (breathing) modes for these nanostruc- 
tures along with those of two higher-order modes. For metal nanoshells, we find that the breathing 
mode frequency is always lower than the one for solid particles of the same size, while the damping 
is higher and increases with reduction of the shell thickness. We identify two regimes that can be 
characterized as weakly damped and overdamped vibrations in the presence of external medium. For 
bimetalUic particles, we find periodic dependence of frequency and damping rate on the shell thick- 
ness with period determined by mode number. For both types of nanostructures, the frequency 
of higher modes is nearly independent of the environment, while the damping rate shows strong 
sensitivity to outside medium. 
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I. INTRODUCTION 

Composite nanostructures comprised of alternating 
metal and dielectric constituents such as, e.g., bimetal- 
lic spherical particles, metal nanoshells,^ nanorice,^, 
and structures with more than two constituents (nano- 
matryoshka)2i, are among the highlights of nanomaterials 
with versatile optical and mechanical properties. Similar 
to solid metal nanoparticles, the absorption and scat- 
tering of light by composite nanostructures are domi- 
nated by the surface plasmon resonance (SPR). These 
structures offer new possibilities for various applications 
due to tunability of their optical properties which can 
be controlled during the manufacturing process by vary- 
ing constituent materials or geometry.**'^ For example, by 
changing shell thickness, the SPR position in nanoshells 
and bimetallic particles can be changed in wide frequency 
range. Recent measurements on single nanoshells^'^ in- 
dicated enhanced sensitivity of their scattering spectra 
to the environment, while the SPR redshift resulted in 
a narrowing of their resonance lineshape as compared 
to solid gold particles.- These unique optical charac- 
teristics of composite structures spurred numerous pro- 
posals for various applications, e.g., in optomechanicS)^ 
sensingjiSiii drug deliver y ^^i^^ and phototherapy^^ 

While the optical response of core-shell nanostruc- 
tures has been extensively studied, much less is known 
about their acoustical properties J^ii^iii In fact, the low 
frequency vibrational modes of composite nanoparticles 
bear a unique signature of their structural and mechan- 
ical characteristics directly reflecting the impact of con- 
finement on the ionic movement j ^^i^^i^^ This is in con- 
trast to the optical response which is mainly determined 
by electronic excitations. Thus the vibrational modes 
constitute an additional source of information that is 
particularly important for hybrid or layered systems with 
structurally distinct constituents, such as, e.g., bimetalhc 
nanoparticles and nanoshells. Recent pump-probe spec- 



troscopy measurements of fundamental (breathing) mode 
frequency and damping in bimetallic particles22i2ii^ and 
gold nanoshella^^ revealed significantly different, from 
solid nanoparticles, acoustical signature that could be 
used as a tool complimentary to optical methods. 

In this paper we study the spectrum of radial vibra- 
tional modes of core-shell nanoparticles, including metal 
nanoshells with dielectric core and bimetallic spheri- 
cal particles, embedded in an environment. Acoustic 
vibrational modes in nanoparticles can be impulsively 
launched by a laser pulse and monitored using the pump- 
probe spectroscopic techniquej ^^i^^'^^i^^'^°i^^ After initial 
stage of rapid expansion, the nanoparticle undergoes ra- 
dial oscillations around its new equilibrium size. The 
periodic change in volume translates into a modulation, 
in time, of the SPR energy and manifests itself in the 
coherent oscillations of differential transmissioni^ Since 
the laser spot is usually much larger than the nanopar- 
ticle size, the initial expansion is homogeneous, so that 
predominantly the fundamental {n — 0) breathing mode, 
corresponding to oscillations of the nanoparticle volume 
as a whole, is excited. In solid nanoparticles, frequencies 
of radial modes are inversely proportional to the nanopar- 
ticle radius, u;„ « 7r(n + 1)clI R^^ where cl is longitu- 
dinal sound velocity. When a nanoparticle is embedded 
in dielectric medium, the oscillations are damped, due 
to the transfer of latice energy to the environment, with 
damping rateS^i^i 7 ex R^^ . For homogeneous intial ex- 
pansion, the higher vibrational modes are much less pro- 
nounced. With the pump-probe technique, the n — 1 
and n = 2 modes were recently detected in Ag spherical 
particles, while the higher radial modes (up to n = 4) 
were observed using Raman scattering technique.^* 

In this paper we present a detailed analysis of en- 
ergy spectrum for the first three {n — 0, 1, 2) lowest ra- 
dial vibrational modes in metal nanoshells and bimetal- 
lic core-shell nanoparticles in different environments. We 
find that the modes eigen-energies exhibit a strong de- 



pendence on the aspect ratio, k = R1/R2, where i?i 
and i?2 are the inner and outer radii, respectively. For 
nanoshells, the fundamental mode energy is considerably 
lower than that for solid particles, while the damping 
rate exhibits a more complicated behavior; at k « 0.4 
it has a minimum, while for larger k, the damping rate 
sharply increases, 7 oc d^^, where d is the shell thickness. 
For thin shells, we find two regimes, weakly damped and 
overdamped, characterized by the nanoshell mass relative 
to displaced mass of the outer medium. These regimes 
are absent for higher modes, n > 1, where both fre- 
quency and damping rate are inversely proportional to 
the shell thickness. At the same time, the impact of envi- 
ronment on higher modes frequencies is negligible, while 
the damping rates are strongly affected by outside dielec- 
tric. For thin shells, the fundamental mode carries most 
of the oscillator strength, implying enhanced oscillations 
amplitude as compared to solid particles. However, the 
higher modes have relatively larger oscillator strengths 
and thus should be better detected. 

For bimetallic Au/Ag nanoparticles, we find that, in a 
wide range of aspect ratios, the fundamental mode fre- 
quency and damping rate are higher than those for ei- 
ther Au or Ag solid particles. For higher modes, both 
frequency and damping exhibit periodic oscillations with 
aspect ratio, in contrast to nanoshells, with period deter- 
mined by the mode number n. Similar to nanoshells, for 
higher modes, the environment significantly affects only 
the damping rate. The oscillator strengths also exhibit 
periodic dependence on k and have maximum values for 
thin shells. 

The paper is organized as follows. In section [TTl equa- 
tions determining the complex eigenvalues for radial vi- 
brational modes of spherical core-shell nanoparticles and 
expressions for their oscillator strength are derived. The 
analysis of radial modes for Au nanoshells with dielec- 
tric core is presented in section IIIII In section IIVI the 
results for bimetallic Au/Ag core-shell nanoparticles are 
discussed. Section IVl concludes the paper. 



II. SPECTRUM OF RADIAL VIBRATIONAL 
MODES OF A CORE-SHELL NANOPARTICLE 

We consider radial normal modes of a spherical core- 
shell particle with core extending up to inner radius Ri 
and shell up to outer radius i?2 embedded in a dielectric 
medium. At zero angular momentum [l — 0), the mo- 
tion of nanoparticle boundaries is described by the radial 
displacement u{r) — ip'{r) (prime stands for derivative), 
where displacement potential ip satisfies the Hclmholtz 
equation, 
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where k^^' = w/c]* is the wave-vector and longitudinal 
sound velocity c}* in ith medium. The boundary condi- 
tions demand continuity of both the displacement u{r) 



and the radial diagonal component of the stress tensor, 
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across the core/shell and shell/medium interfaces. Here 
p and cl,t are the corresponding density and longitudi- 
nal/transverse sound velocities. In the core, shell, and 
medium regions, solutions of Eq. ([T|) are, respectively, of 
the form 
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where is phase shift. Matching u ~ ip' and arr at the 
inner radius i?i and outer radius i?2 yields equations for 
eigenvalues 
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where ^ = k'^'^^R2 = ujR2/cj^ and k = R1/R2 are short- 
hand notations for the normalized eigenfrequencies and 
aspect ratio, respectively. The parameters 
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characterize the metal/dielectric interfaces (i — c,s,m 
stand for core, shell, and outer medium). The parameter 
Xm = 4/3^(1 — Pml l^-s) can be expressed via the ratio of 
shear modules of outside dielectric p,m and metal shell 
Ps, where shear module is defined as p = pc^. Equa- 
tions ([4]) determine complex eigenvalues ^n = ^^ + i^n 
of composite core-shell particle, n = 0, 1, . . . being mode 
numbers, where the real and imaginary parts describe 
mode frequency and damping rate. 

From Eqs. (|4]), the known results for solid particles can 
be easily recovered by setting (p = 0. In the case when 
density of outside medium is much smaller than that of 
particle, the approximate solution for eigenvalues reads 
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where Zi — pic]^ is the acoustic impedance of the ith ma- 
terial. Note that for solid nanoparticles, both frequency 
uj = ^'cl/ R and damping rate 7 = ^' c^jR decrease with 
particle size and that the latter is nearly independent of 
mode number n. 

Let us now turn to relative contributions of different 
eigenmodes in the total acoustical response. We assume 
that the laser spot is larger than the nanoparticle size 



and, therefore, the initial expansion is isotropic. Even 
so, as the nanoparticle vibrates with respect to its new 
equihbrium size, all radial modes are, in principle, ex- 
cited. A precise determination of modes contibutions is 
somewhat complicated due to the fact that, for composite 
heterogeneous nanostructures, these modes are not nec- 
essarily orthonormal. However, since the dominant con- 
tribution comes from the fundamental (breathing) mode, 
the orthonormal set can be constructed perturbatively. 
Namely, the higher modes acquire an admixture of the 
fundamental one, while the latter, in this approximation, 
remains nearly unperturbed. The orthonormal displace- 
ments set has the form 
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characterizes the admixture. The radial displacement 
dr{t) (relative to average radius) can be expanded over 
Vn{r) as 6r{t) — X]n^»^»('')^~*""*> where coefficients 
bn are determined from the initial condition. For ho- 
mogeneous initial expansion, Sr(0) ex r, the oscillator 
strengths, 
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determine relative contributions of the corresponding 
modes in the total acoustical response. In the following 
sections, we present detailed analysis of the vibrational 
modes spectrum for metal nanoshells and bimetallic par- 
ticles. 



III. VIBRATIONAL MODES OF METAL 
NANOSHELL 

The approach developed in previous section describes 
vibrational eigenmodes of a nanoparticle in an equi- 
librium state. However, acoustical modes are usually 
launched by ultrashort laser pulse that results in rapid 
initial expansion of the particle followed by oscillations 
around its new equilibrium position. In the case of a 
metal nanoshell with a dielectric core, the latter plays 
negligible role and the eigenmodes are determined mainly 
by the metal shell and outside medium. Indeed, the di- 
electric core is not directly affected by the laser pulse, 
but experiences thermal expansion as a result of heat 



transfer from the metal shell. At the same time, this ex- 
pansion is much weaker than that of the metal, so that 
when new equilibrium size is established, the core is al- 
most completely disengaged from the shell. This should 
be contrasted to bimetallic particles where the core does 
remain engaged after the expansion and thus contributes 
to the acoustical vibration spectrum (see the next sec- 
tion). 

Consider first the fundamental {n — 0) mode. For the 
ideal case of a nanoshell in vacuum, frequency is obtained 
from Eq. ^ with stress-free inner boundary, i.e., by set- 
ting rjm = rjc = and Xc — Xm = 4/3^. In the thin shell 
limit, 1 ~ K — d/i?2 ^ 1, we then recover the known 
result^ 



^0 = 2/3,73^4^ 
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For a nanoshell in a dielectric medium, the eigenvalues 
are complex due to the energy exchange with the envi- 
ronment. For thin nanoshells, 1 — k = d/R2 <^ 1, Eqs. 
(|4|) (with rjc — 0) reduces to 
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In the typical case when the metal shell density is much 
higher than that of the surrounding medium, i.e., ?],„ = 
p(™)/p(s) < 1, Eq. dni) further simplifies to 
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where x = S,/S,o and we used x 
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For thin shells, we now can identify two regimes gov 
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where Mg is the metal shell mass, and Mm is the mass 
of outer medium displaced by the core-shell particle. For 
a "heavy shell", Ms ^ M,„, the complex eigenvalue is 
given by 
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where ^0 is the eigenvalue for a nanoshell in vacuum and 
A — amilm/S,o{l — k). In a good approximation, the real 
part is simply ^' ~ ^qj and is thus independent of the 
medium or aspect ratio. In contrast, the imaginary part, 
however small (^" <C (J), is only non-zero in the presence 
of outside medium and thus depends on both. Putting 
all together, we obtain in the "heavy shell" regime 



,(0) 



.(0) 



2c^ 

i?2 



V^^Ws 



277^/3,^(3 -4/3g) 
d a^ -f 4/32(3 -4/32)- 



(17) 



TABLE I: Longitudinal velocity cl (ni/s), transverse velocity 
CT(ni/s), density p(kg/ni'') acoustic impedance Z(10'^kg/m'^s) 
and shear module /i(10^"kg/ms^) of bulk gold, silver and sev- 
eral dielectrics. 
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Au 


3240 


1200 
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6.38 
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3.80 


2.87 


Glass 


4610 


2590 


3810 


1.76 


2.56 


Fused silica 


5970 


3765 
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1.31 


3.12 


Water 


1490 





1000 


0.15 






Here the damping rate is determined by the shell thick- 
ness rather than by the overall size. In the opposite case 
of a "light she ll" , Ms <C Af™, the eigenvalue is given by 
^ ~ 2am/3m(v^l -/3^ + i/3m), yielding 
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Note that, in this case, if the transverse velocity of out- 
side medium is zero (e.g., in water), both frequency and 
damping rate vanish. 

Below we present the results of numerical calculations 
of vibrational modes spectrum for Au nanoshells in three 
dielectric media: water, fused silica and glass (BaO- 
P2O5). The material parameters are listed in Table [H 
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FIG. 2: Normalized frequency (a) and damping rate (b) of 
n — 1 radial mode in Au nanoshell. Dashed line in cartoon 
indicates zero-displacement surface. 
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FIG. 1: Normalized frequency (a) and damping rate (b) of 
fundamental (n — 0) radial mode in Au nanoshell vs. aspect 
ratio in various environments. The cartoon illustrates the 
motion of shell boundaries. 



In Fig. [T] we show the energy and damping rate for fun- 
damental {n — 0) breathing mode versus aspect ratio 
i?i/i?2- In entire range of aspect ratios, the nanoshell 
frequencies are lower than those of solid particles of the 
same overall size. Both "heavy" and "light" regimes can 
be clearly distinguished. In the former regime, for as- 
pect ratios k < 0.5, the eigenfrcquencies are only weakly 
dependent on media and follow each other, according to 
Eq. (fT7)) . For aspect ratios k > 0.5, the eigenfrequencies 
start deviating from each other due to the effect of en- 
vironment, and for very thin shell (k > 0.9) the "light 
shell" behavior is observed, in agreement with Eq. (llSp . 
In contrast, the damping rate strongly depends on the en- 
vironment for all K, however, the two regimes are clearly 
seen. For thick shells, the damping is mostly determined 
by longitudinal sound velocities of outside media [see 
Eq. ([T7|) ]. while for thin nanoshells it is dominated by 
the transverse ones [see Eq. (IT5] so that the curves for 
glass and fused silica cross each other. Importantly, in all 
cases, the damping rate exhibits a minimum at k « 0.4, 
while decreasing shell thickness it behaves as 7 ex d~^. 
Note that for thin nanoshell in water, both frequency 
and damping rate vanish due to the absence of trans- 
verse sound (cy"*'^'' = 0). 

Situation is drastically different for higher modes, 
shown in Fig. [2] and [3] for n = 1 and n = 2, respec- 
tively. Both frequency and damping rate increase with 
aspect ratio and behave as d~^ for thin shells. The reason 
is that higher modes correspond to oscillations around 
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FIG. 3: Same as Fig. [H but forn = 2 mode. 



a zero-displacement surface inside the shell rather than 
modulation of its overall size. In this case, the cigcnfre- 
quencies are nearly independent of environment in the 
entire range of aspect ratios, indicating the presence of 
only one ("heavy" shell) regime. At the same time, the 
damping rate exhibits a weak minimum that shifts to- 
wards smaller k as n increases (compare Figs. [Hand [3|). 

In Fig. m we show calculated oscillator strengths for 
n — 0,1,2 modes versus aspect ratio for a nanoshell in 
vacuum. In solid nanoparticles, the oscillator strength for 



n-mode can be approximated by C„ 



[sol) 



3V10/[7r(n 



1)]^. However these values slightly larger than those ob- 
tained numerically (k = in Fig. 3]). This discrepancy 
can be attributed to the corrections due to the second 
term in Eq. ([6]) (third term vanishes in vacuum) that re- 
duce the oscillator strength of, e.g., n = 1 mode from 0.25 
to 0.2. In contrast, for nanoshells, C; vanishes in the thin 
shell limit for all higher modes I ^ 0. For 5 = d/ R2 ^ 1, 
the behavior of C„ can be approximated as 
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In the d — limit, Cq reaches it maximal value, Cq = 
1, i.e., the fundamental mode carries the entire oscil- 
lator strength. As a result, for thin nanoshells, the 
fundamental mode is more pronounced as compared to 
solid particles.— At the same time, for thick nanoshells, 
R1/R2 « 0.3, the oscillator strength of n = 1 mode is 
considerably higher than in solid particle (see Fig. [1]). 
This enhancement is less pronounced for higher modes. 
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FIG. 4: Oscillator strengths for n = 0, 1 (a) and n — 2 (h) 
radial modes of a nanoshell in vacuum vs. aspect ratio is 
shown together with analytical fit. 



IV. VIBRATIONAL MODES OF BIMETALLIC 
NANOPARTICLE 

Consider now bimetallic core-shell nanoparticles em- 
bedded in outside dielectric. As discussed in the previ- 
ous section, the acoustical modes are determined by vi- 
brations of both core and shell parts. Below we present 
results for eigenfrequencies and damping rates of Au/Ag 
core-shell spherical particles in different environments 
based on full numerical solution of Eqs. ^. Since com- 
plex eigenvalues ^„ are normalized in terms of longitudi- 
nal sound velocity in the shell region, in order to obtain 
eigenvalues for solid Au sphere at k = 1 one simply needs 
to rescale ^„ by c\^ jc^h ■ 

In Fig. [51 we show the frequency and damping rate 
vs. aspect ratio in three different media. In contrast to 
nanoshells, here the frequency is higher than in either Ag 
and Au solid particles in a wide range of k. Another dis- 
tinction is a dependence of the breathing mode frequency 
on the environment in entire range of k. The sensitivity 
to the environment is much more pronounced for damp- 
ing rate that shows a weak maximum for aspect ratios 
K w 0.5. In the course of transition from solid Ag to 
solid Au nanoparticles, the magnitude of damping drops 
roughly by factor of two reflecting the relation between 
acoustic impedances of Ag and Au (see Table [J). 

In Figs. [6] and [3 we show the result of calculations for 
n = 1 and n = 2 modes, respectively. The main fea- 
ture distinguishing bimetallic particles from nanoshells 
is periodic dependence of both frequency and damping 
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FIG. 5: Frequency (a) and damping rate (b) of n = radial 
mode in Au/Ag bimetallic particle, normalized to the shell 
longitudinal sound velocity, in various environments. 



rates on aspect ratio. The number of oscillations coin- 
cides with the mode number n. Similar to nanoshells, the 
frequencies of higher modes only weakly depend on the 
environment, while the damping rate shows the opposite 
trend. Furthermore, the oscillation amplitude of damp- 
ing rate increases with aspect ratio reaching about 50%. 
The origin of oscillations is the interference of cavity and 
shell modes with change of relative core and shell sizes. 

Oscillator strengths of n = 0, 1, 2 modes of Au/Ag par- 
ticle in vacuum, C„ are plotted in Fig. [5] vs. aspect ratio. 
As in solid particles and nanoshells, here the fundamental 
mode carries most of the oscillator strength, but slightly 
oscillates with aspect ratio. The oscillations are consid- 
erably stronger for Ci and C2, reaching about 50% in 
peak-to-peak amplitude. Note that the higher modes, 
plotted here, are the orthogonal states with amplitude 
Vn (see Eq. [S]) that carry an admixture from the funda- 
mental mode, rather than the pure n = 1 and n = 2 
states. The dependence of admixture coefficients Vn on k, 
is highly non-monotonic, reaching maximal value about 
0.3 for n — 1 mode and 0.2 for n = 2 mode (see inset). 
This admixture is reducing the strength of n = 1 mode, 
and enhancing that of n = 2 mode, resulting in about the 
same oscillator strength for both modes for aspect ratio 
close to 0.4. Note, finally that, in contrast to nanoshells, 
higher modes should be more detectable for higher aspect 
ratios (k k, 0.8 for Au/Ag particle). 



FIG. 6: Same as Fig. (5] but with n — 1. 
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FIG. 7: Same as Fig. [S] but with n = 2. 



V. CONCLUSIONS 

In summary, we studied the spectrum of radial vibra- 
tional modes in bimetallic core-shell nanoparticles and 
metal nanoshells in an environment. We found that the 
acoustical signature of these nanostructures is extremely 
sensitive to their composition, geometry and environ- 
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FIG. 8: Oscillator strengths are shown for n — (a.) and 
n — 1,2 (b) radial modes of Au/Ag particle in vacuum vs. 
Ri/R2- Inset: admixture coefficients for n — 1,2 modes. 



ment, and differs significantly from tfiat of solid nanopar- 
ticles. In particular, for bimetallic nanoparticles, fre- 
quency and damping rate are periodic functions of the as- 
pect ratio, R1/R2, with period determined by the mode 
number, while in nanoshells this dependence is mono- 
tonic. This difference between nanoshels and bimetallic 
particles can be traced to the fact that, in the latter, the 
core is fully engaged in vibrations, leading to the inter- 
ference between cavity and shell modes. Another distinc- 
tion is much higher damping rate in thin nanoshells in 
comparison with solid or bimetallic particles. Thus, the 
acoustical probe of composite nanostructures can provide 
a sensitive tool for determination of their composition 
and geometry . 
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